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ABSTRACT 

We study the excitation of nonaxisymmetric modes in the post-merger phase of binary com- 
pact object mergers and the associated gravitational wave emission. Our analysis is based 
on general-relativistic simulations, in the spatial conformal flatness approximation, using 
smoothed-particle-hydrodynamics for the evolution of matter, and we use a set of equal and 
unequal mass models, described by two nonzero-temperature hadronic equations of state and 
by one strange star equation of state. Through Fourier transforms of the evolution of matter 
variables, we can identify a number of oscillation modes, as well as several nonlinear com- 
ponents (combination frequencies). We focus on the dominant m = 2 mode, which forms a 
triplet with two nonlinear components that are the result of coupling to the quasiradial mode. 
A corresponding triplet of frequencies is identified in the gravitational wave spectrum, when 
the individual masses of the compact objects are in the most likely range of 1.2 to 1.35 Mq. 
We can thus associate, through direct analysis of the dynamics of the fluid, a specific fre- 
quency peak in the gravitational wave spectrum with the nonlinear component resulting from 
the difference between the m = 2 mode and the quasi-radial mode. Once such observations 
become available, both the m = 2 and quasiradial mode frequencies could be extracted, al- 
lowing for the application of gravitational-wave asteroseismology to the post-merger remnant 
and leading to tight constraints on the equation of state of high-density matter. 
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I 1 INTRODUCTION 

Mergers of binary compact objects are prime sources for second- 
and third-g eneration interfer o metric gravitational w a ve (GW) de- 
tector s jAbbott et all |2009| ; lAcernese etail l2006t lAbadie et all 
l20ld) . The expected GW signals from such events are estimated 
throu gh general-relativistic hydrodynamical simulations (see lDuezl 
|2010| for a review). The outcome of these simulations depends on 
the binary parameters and on the equation of state (EOS) of high- 
densi t y matter. The latter is rather uncertain (lLattimer & Prakasn 
l2007l : lsteineret al.ll2010h and currently it is unclear whether, for 
a given binary mass, the merger would lead to a hypermassive 
compact object or to prompt collapse to a black hole. If a hy- 
permassive compact object forms, it will not be axisymmetric, 
but it will show transient nonaxisymmetric deformations, such as 
a bar-like shape, spiral arms, a double core structure and quasi- 
radial and nonaxisymmetric oscillations of the matter. Strong non- 
axisymmetric features should be distinguishable in the GW signal 
and could be used for characterizing the hypermassive compact ob- 
ject fehuge et al.ll99i:l0echslin et al.l2002l;IShibata & Urvul2002l; 
Shibata et alj|2005l; IShibata & Taniguchill2006l: lOechslin & Jankal 
2007l:lBaiottiet"al. 2008; Kiuchi et alj2009; Bauswein et alj2010h. 



Here, we analyze the formed hypermassive compact object as 
an isolated gravitating fluid, studying its oscillation modes. Fourier 



transforms of the evolved variables reveal that the fluid is oscillat- 
ing in a number of modes that have discrete frequencies throughout 
the star. Furthermore, we also identify several nonlinear compo- 
nents, sums and differences of discrete oscillation modes. The os- 
cillations identified in the fluid are in direct correspondence with 
peaks in the GW spectrum, as obtained through the quadrupole 
formula. We focus on the main quadrupole (m = 2) oscillation 
mode of the fluid, which appears as a triplet, the side bands be- 
ing due to the nonlinear coupling to the fundamental quasiradial 
(m = 0) mode. The lowest-frequency side band, the difference be- 
tween the m — 2 and m — frequencies, coincides with a peak in 
the GW spectrum that characterizes the merger phase in all mod- 
els, in which the mass of both compact objects are in the range of 
1.2 to 1.35M0. We thus propose that in the event of detection, one 
could extract both the m — and m — 2 mode frequencies of the 
merged object, which could lead to tight constraints on the EOS 
of high-density matter. In effect, such a detection would enable an 
analysis with GW asteroseismology of the remnants of binary com- 
pact object mergers. 

The simulation code is based on general relativistic smoothed- 
particle-hydrodynamics (SPH) and on a spatially- conformally-flat 
spacet ime approximation, as described in detail in lBauswein et al.l 
(2010) and references therein. We simulate the merger of the com- 
pact objects following the evolution from the late inspiral phase 
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through the merging and the formation of the hypermassive rem- 
nant, until a quasi-stationary state is reached. 

The paper is structured as follows: In Sec. |2]we discuss the 
initial data used in our similations. In Sec.[3]we outline the numer- 
ical methods used in the simulations and in the mode extraction 
procedure. Section [4]presents our results in detail, while in Sec. [5] 
we compare them to previous results in the literature. We conclude 
with Sec.|6] 



2 INITIAL DATA 

We cons ider two differen t models of hadronic EOSs, refer red to 
as Shen dShen et alj[l998l) and LS dLattimer & Swestvl[l99lt) and 
a strange quark matter EOS, referred to as MIT60. Since dur- 
ing merging temperatures of more than several ten MeV can be 
reached, al l three EOSs take int o account non-zero temperature ef- 
fects (see lBauswein et al]|2010l for a discussion of the importance 
of temperature effects in the merger context). 

The Shen EOS was derived within a relativistic mean-field 
theory and assumes an i ncompressibilty m odulus for nuclear mat- 
ter of K = 280 MeV dShen et al.ll 19981) . Solving the relativistic 
equations of hydrostatic equilibrium, nonrotating neutron stars de- 
scribed by the Shen EOS have radii of about 15 km in the mass 
range of 0.5 to 1.7 Mq. The maximum mass of nonrotating stars is 
2.2 Mq. For the LS EOS a liquid-drop model with K = 180 MeV 
was adopted dLattimer & S wests! 199 lh . which yields neutron stars 
more compact in comparison to the Shen EOS . Typical radii are 
of the order of 12 km for masses of 0.5 to 1.6 Mq. The LS EOS 
supports nonrotating objects with masses slightly above 1.8 Mq. 

The MIT60 EOS follows fr om the strange matter hypothe- 
sis dBodmeJ fl97ll ; IWittenl 1 1984) . i.e. that 3-flavor quark matter 
of up, down and strange quarks is more stable than ordinary nu- 
clear matter. In such a case, co mpact stars would be strange stars 
rather than neutron stars (see e.g. lGlendennin g 1996; Haensel et al.l 
2007), which has not been strictly ruled out theoretically or ob- 
servati onally at this point . For the MIT60 EOS the MIT Bag 
model dFarhi& Jaffe1ll984) was employ ed, with a bag constant of 
60MeV/fm 3 (see Bauswei n et alJuOlol for details). Strange star 
models are generally more compact than neutron star models of 
the same mass. Furthermore, the mass-radius relation of strange 
stars does not show the typical inverse relation of neutron stars. 
The mass-radius relatio n for all three EOSs inc luded in our study 
can be found in Fig. 2 of lBauswein et alj J2010h - 

We note that the recent discovery of a 2Mq pulsar 
dDemorest et al.l2010l) practically rules out EOSs which yield max- 
imum masses of nonrotating compact stars below this limit. The 
maximum masses of nonrotating stars for the LS and MIT60 
EOSs employed here fall short of this requirement, but not 
dramatically. Though to date several finite-temperature EOSs 
iTypel et al'j|201Ct iHempel & Schaffner-Bielichll201ol; Ishen et all 



a lypet et 
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2011; Shen et i 



ipei & ^cnatlner-iSielicn ZU1U; J>nen et al. 
al. 2011) have been published, at the time 



of writing only those introduced above have been available to us 
and successfully incorporated in the simulation code. Thus, leav- 
ing out the LS and MIT60 EOSs would leave us with a single EOS, 
which would not test the sensitivity of our results to the choice of 
EOS. We thus include the LS and MIT60 EOSs in our study for the 
sole purpose of estimating the sensitivity of our conclusions to the 
EOS employed, keeping in mind their disadvantage. Nevertheless, 
the individual stars that make up the binaries have masses signifi- 
cantly lower than the maximum mass limit, and the hypermassive 
object formed as a result of the merger is supported against col- 



lapse by strong differential rotation. Because of this, the outcome 
of the simulations with the LS and MIT60 EOSs is still useful for 
qualitative comparisons to the cases where the Shen EOS is used. 

The masses of compact stars in binarie s cluster at 
about 1.35 Mq dThorsett & Chakrabartvl Il999t IZhang et al.l 
which is also predi cted by population synthesis studies 
(e.g. iBelczvnski et al. 2008 . Therefore, we focus in our analysis 
on systems with two neutron stars with a gravitational mass of 
1.35 Mq. In order to investigate unequal-mass mergers, we also 
consider configurations with a 1.2 Mq neutron star and a more 
massive companion of 1.35 Mq. For the MIT60 EOS, we include a 
configuration with two low-mass stars of 1.1 Mq each, an unequal- 
mass binary with 1.2 Mq and 1.35 Mq, and a configuration with 
two stars of 1.35 Mq each. In the latter case, however, the hyper- 
massive object formed after merging only survives for a few dy- 
namical timescales before collapsing to a black hole, and we do not 
discuss this model in our study. Table [T]lists all models considered 
further, where, for example, Shen 12135 refers to the simulation 
with a 1.2 Mq and a 1.35 Mq star employing the Shen EOS. 



3 NUMERICAL METHODS 
3.1 Simulations 

The hydrodynamical simulations of the binary merger and the post- 
merger remnant are performed with a three-dimensional general- 
relativistic SPH code. The Einstein field equations are solved 
within the condition of spatial conformal flatness, which requires 
the additional impleme ntation of a GW backreactio n scheme. For 
details of the code see Oechslin et a I] d2002ll2007h . The calcula- 
tions start from a quasi-equilibrium orbit about three revolutions 
before the actual merger. For the neutron star models about 500,000 
SPH particles are used, while the flat density profile of strange stars 
allows for a lower SPH particle resolution of about 130,000 parti- 
cles. Simulations with a somewhat higher number of particles do 
not significantly affect the main features of our results. The simu- 
lations are carried out until a stable remnant in approximate rota- 
tional equilibrium has formed and oscillates for tens of dynamical 
timescales, or until the delayed collapse to a black hole occurs. 

The dynamics of the models used in th is study have been ex- 
tensively discussed in Bauswein et al. (2010), and the reader is re- 
ferred to this publication for details. Here we only give a brief 
outline: While orbiting around the common center of mass, the 
stars approach each other increasingly faster due to angular mo- 
mentum and energy losses by the GW emission. Prior to merging, 
tidal forces start to deform the stars. The strength of this effect de- 
pends on the EOS and on the masses of the binary components. 
The deformation is pronounced for the Shen EOS, but relatively 
small for the strange star models. Finally, the stars merge into a ro- 
tating double-core structure, where the two dense cores appear to 
bounce against each other a few times to ultimately form a differ- 
entially rotating hypermassive object with a single core. In the case 
of unequal masses, such as our models Shen 12135 and LS 12135, 
the less massive binary component is tidally stretched and wrapped 
around the more massive one. Moreover, an extended spiral arm 
develops, feeding a dilute halo around the differentially rotating 
central object. In the case of the unequal mass strange star model 
MIT60 12135, the outcome is more similar to the equal mass case, 
because of the reduced tidal effects. 

The equal-mass neutron star mergers also form thick disks 
around the remnants, because matter is shed off from the whole 
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surface of the hypermassive objects. In contrast, the remnants of 
merging strange stars terminate at a sharp surface, as did the initial 
stars. Matter is stripped off only later in the evolution when angular 
momentum transfer leads to the formation of two thin spiral arms, 
which then form a fragmented thin disk around the central object. 

The hypermassive remnants of the merger are supported 
against gravitational collapse by differential rotation and thermal 
pressure. However, ongoing angular momentum transfer and loss 
by mass shedding finally leads to gravitational collapse and black 
hole formation. The delayed collapse can take up to hundeds of 
milliseconds for low-mass models, but our simulations were termi- 
nated at about 10 to 20 ms after merging. 

The GW emission of our models is analyzed by means of the 
quadrupole formula. In the low-frequency regime the spectra are 
not reliable because we simulate only a small fraction of the inspiral 
phase and thus miss most of the low-frequency part of the signal. 

3.2 Mode Extraction 

The original simulation data in different two-dimensional planes 
are first mapped onto a Cartesian grid of size 71^X71^=51x51. 
We then perform Fourier analysis of evolved variables over the 
whole grid (but focusing mainly on the equatorial plane) and iden- 
tify the predominantly excited modes, which appear as discrete 
frequencies in the Fourier spectrum. At selected mode frequen- 
cies, we proceed with extracting the s hape of the mode eigen - 
functions, using a method described in Stergioulas etal.1 J2004) . 
The two-dimensional shape of the eigenfunction in the equatorial 
plane correlates directly with the amplitude of the Fourier trans- 
form at a given mode frequency. At nodal lines, care must be 
taken to switc h the sign of the eigenf unction. This method was 
first applied in lstergioulas et all (2004) to the case of axisymmet- 
ric modes, where two-dimensional eigenfunctions were extracted 
in a plane that passes through the rotational axis. Here, we find that 
the same method can be used for extracting the eigenfunctions of 
non-axisymmetric modes in the equatorial plane. We note that for 
all modes that were identified, the mode frequency was discrete 
throughout the star and identical in all hydrodynamical variables. 

To our knowledge, this is the first direct demonstration that 
once a hypermassive object forms after a binary merger, it behaves 
as an isolated, self-gravitating system and its dynamics can be de- 
scribed as a superposition of different oscillation modes (in a differ- 
entially rotating background). Departures from a background equi- 
librium state can be described (up to a certain degree) by nonlinear 
features of the oscillations. We note that the background is varying 
in time, so that individual frequency peaks are broadened. Never- 
theless, the frequency peaks remain sufficiently narrow to identify 
individual modes. 
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Figure 1. Projection in the equatorial plane of the extracted m = 0, 2 and 
4 mode eigenfuctions (from top to bottom) for the oscillations in pressure 
of model Shen 135135. The two axes count individual grid points of the 
Cartesian grid used for the mode analysis. The color scale only has a relative 
meaning. 
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Figure 2. Projection in the equatorial plane of the extracted m = 1 mode 
eigenfunction, for the oscillations in pressure of model MIT60 12135. 



stars. The two nodal lines are, as expected, nearly perpendicular to 
each other. The m = 4 mode eigenfunction shows a characteris- 
tic octupole pattern. The alternating regions are of nearly the same 
size, but due to the very low amplitude of excitation and the finite 
duration of the time-series, the Fourier transform also picks up a 
non-zero contribution at the center of the star, which means that for 
this mode we are approaching the accuracy of our method. Higher- 
order modes should also be excited (and there are indeed a large 
number of peaks in the Fourier transform at high frequencies) but 
their amplitude is so small that any extracted eigenfunction will be 
dominated by numerical errors. 

Fig. [2] shows the eigenfunction of the m = 1 mode excited 
during the merger of the uneqal-mass model MIT60 12135. A nodal 
line that cuts through the center of the star is characteristic of this 
mode. Notice that conservation of momentum prohibits the exis- 
tence of a fundamental m = 1 mode, as this would induce a linear 
motion of the center of mass. However, higher-order overtones are 
allowed, so the extracted mode should be identified as the (lowest- 
frequency) first overtone of the m — 1 modes. This is supported by 
the additional, nearly circular node line, which allows for momen- 
tum conservation during the oscillation. 



4 RESULTS 

4.1 Two-dimensional eigenfunctions 

Non-isentropic, differentially rotating, compact objects can oscil- 
late in a large number of different quasinormal modes, (/-modes, p- 
modes, g-modes and inertial modes, distinguished by three differ- 
ent mode numb e rs, I, m and n, the latter one being the radial order, 
see I Stergioulas! |2003) for a review)). Here, we only characterize 
the structure of the extracted modes in the equatorial plane, in terms 
of the mode number m. A complete characterization and identifica- 
tion with specific quasinormal modes would require a much more 
detailed analysis that is out of the scope of the current work. In par- 
ticular, neither a linear, nor a nonlinear complete analysis of pos- 
sible oscillation modes has been done for remnants of mergers of 
compact objects, to date, with which we could directly compare 
our results. As such, we will talk about "modes" and "eigenfunc- 
tions" even though we have not strictly shown a correspondence to 
specific linear quasinormal modes. Nevertheless, for the purpose of 
explaining the main observable features of GWs from nonaxisym- 
metric oscillations in binary compact object mergers, our present 
approach suffices. 

In the Fourier transform of hydrodynamical variables, the 
dominant mode is always the m = 2 mode. For equal-mass merg- 
ers, the even mode numbers are predominantly excited, while for 
unequal-mass mergers, the odd mode numbers are also consider- 
ably excited (mainly the m = 1 mode). The modes can be identi- 
fied by inspection of the extracted two-dimensional eigenfunction 
in the equatorial plane. Representative cases are shown in Figs.Q] 
and [2] which display the eigenfunctions (in pressure oscillations) 
of the m = 0, 2 and 4 modes for the equal-mass model Shen 
135135 and the m = 1 mode for the unequal-mass model MIT60 
12135. The m = (quasiradial) eigenfunction is dominated by 
a spherically symmetric contribution, while an additional, small 
(rotationally-induced) quadrupole distortion is also present. As ex- 
pected for this type of eigenfunction, there is a nearly circular nodal 
line in the pressure perturbations, characteristic of the fundamental 
radial mode. 

The eigenfunction of the m — 2 mode shows a very sharp 
quadrupole pattern, induced by the merger of the two individual 



4.2 Mode frequencies 

Table [T] lists the extracted mode frequencies (in the inertial frame, 
defined by coordinate time t) for the various models. We note 
that the frequencies of the quasiradial (m = 0) mode differ sig- 
nificantly among the three chosen EOSs. In the selected mass- 
range, this frequency is lowest for the Shen EOS (about 0.5kHz). 
It doubles to about 1kHz for the LS EOS and it reaches up to 
about 1.5kHz for the MIT60 EOS. This mode frequency is affected 
mainly by two factors: the compactness of the star and the distance 
of the model from the region of axisymmetric instability to col- 
lapse. At the boundary of this region, the frequency of the quasi- 
radial mode goes through zero. On the other hand, higher compact- 
ness leads to higher mode-frequency. But, because of the axisym- 
metric instability, the quasiradial mode frequency can be small for 
any EOS, as long as the model is near the instability threshold. This 
latter property explains the fact that, for each EOS, the models with 
larger total mass have smaller quasiradial mode frequency. 

The frequency of the m = 2 mode depends mainly on the 
compactness of the star and for the models in Table[T]it is roughly 
2 kHz for the stiffer Shen EOS and roughly 3kHz for the softer 
LS EOS and for the MIT60 EOS. As expected, for each EOS, it 
is higher for the more massive (hence more compact) model. In 
all unequal mass cases, we were able to extract the m = 1 mode 
frequency, which is induced by the unequal distribution of mass just 
prior to merging. The additional m = 3 mode was extracted for the 
Shen and LS EOSs, while in three models we could also identify 
the m = 4 mode. 

A remarkable property of the nonaxisymmetric modes m = 
1...4 is that the extracted frequencies are nearly integer multiples 
of the frequency of the m = 1 mode. This property holds true with 
good accuracy for the equal-mass cases, while for the unequal-mass 
mergers it is less accurate, but still roughly true. 



4.3 A representative model 

The evolution of the GW amplitude h+ expected from the late bi- 
nary inspiral, merger and post-merger phases for the unequal-mass 
model Shen 12135 is shown in the top panel of Fig. [3] The merger 
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Table 1. Extracted mode frequencies. 



Model 


/m=0 


/m = l 


fm=2 


fm — 3 


fm=4 




(kHz) 


(kHz) 


(kHz) 


(kHz) 


(kHz) 


Shen 12135 


0.50 


1.07 


2.14 


3.22 


4.26 


Shen 135135 


0.46 




2.24 




4.11 


LS 12135 


1.10 


1.55 


3.12 


4.66 




LS 135135 


0.98 




3.30 






MIT60 1111 


1.56 




2.92 




5.94 


MIT60 12135 


1.24 


1.72 


3.26 







takes place at around t = 17ms from the start of the simulation. 
A characteristic sudden change in the frequency and phase of the 
waveform takes place at the onset of merging. After a few highly 
nonlinear oscillations, the waveform settles into a more regular 
pattern. The late-time damping of the waveform is the result of a 
combination of shock-induced dissipation of oscillations and nu- 
merical damping. Gravitational-wave damping is effective only on 
larger timescales than shown here. With less numerical damping, 
the gravitational wave amplitude could remain strong for a much 
longer time, leading to an enhance d signal (for a recent long-term 
evolution see Rezzoll aet al.l d2010h ). 

Irrespective of the true damping timescale, we are inter- 
ested mainly in the frequency spectrum, which is shown in the 
middle panel of Fig. [3] We plot the scaled power spectral den- 
sity, /i+(/)y / J (black curve), where h+(f) is the Fourier trans- 
form of h+(t), which is directly comparable with the antici- 
pated sensitivity for Advanced LIGO pre s ently being installed 
Marry & the LIGO Scientific Collaboration! (I2010|)) and t he pro - 
jected Einstein Telescope (ET) detectors dPunturo et al.l d2010h : 
Ifflld et alj j2O10h ) (dotted red and blue curves, res pectively). Note 
that th e sensitivity of the advanced Virgo detector dAcernese et al.l 
(2006)) is planned to be similar to the one of the Advanced LIGO 
detector. Therefore, we include the latter only as an example for the 
class of detectors going into operation within the next years. The 
distance to the source is assumed to be 100 Mpc. We also show 
the same quantity during the pre-merger phase only (red curve) and 
during the post-merger phase only (green curve). The first peak at 
roughly 900Hz is artificial, because our simulations only start at a 
few orbits before merging, so most of the actual inspiral phase is 
missing. The largest peak (denoted by /2) at about 2kHz is clearly 
produced exclusively in the post-merger remnant. In addition to this 
peak, which is very pronounced and clearly detectable by ET and 
with good prospects for detection by Advanced LIGO, there are 
two more characteristic peaks, denoted by /_ and /+, which are 
also produced in the post-merger phase and have good prospects 
for detection by ET (more so for /_ than for /+). These two peaks 
are both above the ET noise curve and above the contribution of the 
pre-merger (inspiral) waveform to the FFT of the signal. There are 
additional peaks above the ET noise curve in the range of 1kHz to 
1.5kHz, which, however, will be superseded by the inspiral signal, 
unless the post-merger signal is analyzed separately from the inspi- 
ral part. Furthermore, the high-frequency peaks that are below the 
ET noise curve will not concern us. 

The three frequencies /_ , fe and /+ form a triplet, which we 
will attempt to interpret in terms of the oscillation modes of the 
post-merger remnant. The bottom panel of Fig. [3]shows the ampli- 
tude of the Fourier transform of the evolution of the pressure along 
a fixed coordinate line in the equatorial plane which passes through 
the center of the coordinate system (this is an integrated ampli- 
tude, taking into account the contributions along the whole coordi- 



nate line). The Fourier amplitude is composed of several very dis- 
tinct and narrow peaks, which correspond to the discrete m = 0, 
1, 2, 3 and 4 mode frequencies. In addition, a number of smaller 
peaks are seen. The latter are nonlinear components, combination 
frequencies of th e main oscil l ation m odes, i.e. linear sums a nd dif- 
ferencefl (see IZanotti et all fc005h ; IPassamonti et alj fc007l) and 
references therein). The difference of the frequencies of the m = 2 
and m = mode is denoted as "2-0", their sum as "2+0" (and sim- 
ilarly for other components). A two-dimensional plot of the Fourier 
amplitude of the pressure evolution at the "2-0" frequency reveals 
that indeed it is a combination of a radial and a quadrupolar pattern. 

By comparison of the middle and bottom panels of Fig. [3] 
there is an obvious near coincidence of the frequency peak fa in 
GWs with the frequency of the m = 2 mode of the post-merger 
remnant. Furthermore, the other two peaks of the GW triplet, /_ 
and /+, nearly coincide with the "2-0" and "2+0" combinations of 
the pressure oscillations. As we will see, this is a generic behaviour 
for all models (except for the very low mass MIT60 1111 model). 
It is thus tempting to attribute the /_ peak in the GW spectrum 
to a nonlinear interaction between the quadrupole and quasiradial 
modes. Essentially, a double core structure that appears and dis- 
appears several times in the early post-merger phase could be the 
result of (or could be described by) this nonlinear interaction. To 
firmly establish this goes beyond the scope of the present work, but 
already the coincidence of the /_ and "2-0" frequencies has prac- 
tical consequences: it allows the determination of both the m = 2 
and m = oscillation frequencies of the post-merger remnant. 
In combination with the additional knowledge of the characteris- 
tics of the compact obj ects from the detection of the inspiral signal 
(see lRead et alj |2009i) and references therein), the determination 
of at least two independent oscillation frequencies should allow for 
the application of the theory of GW asteroseismology, placing very 
stringent constraints on several properties of the post-merger rem- 
nant and consequently on the EOS of high-density matter (for an 
applicati on of GW asteroseismolo gy to isolated, rotating neutron 
stars, see lGaertig"& KokkotaJ j201lr) ). 

4.4 A survey of different models 

Having examined in detail model Shen 12135 as a representative 
case, it is instructive to compare the main results to the equal-mass 
case Shen 135135 (Fig.[4]l. The GW spectrum (middle panel of Fig. 
|4]l is very similar to the representative case, except for small differ- 
ences in the frequencies of the main peaks. In the Fourier transform 
of the pressure evolution (bottom panel of Fig. |4), the odd modes 
m — 1 and 3 are no longer as strong as in the unequal-mass case, 
and the evolution is mainly determined by the m = 2 mode, while 
the m — and m — 4 modes are also present and so are nonlinear 
components. Again, there is a coincidence between the frequency 
of the /_ peak in the GW spectrum and the "2-0" combination fre- 
quency in the pressure evolution. 

The unequal-mass model LS 12135 shows similar qualitative 
properties as the corresponding Shen 12135 model, except that 
all mode frequencies are considerably higher. In addition, a larger 
number of nonlinear components can be clearly identified (bottom 
panel of Fig.|5j. The /_ vs. "2-0" coincidence remains. The cor- 
responding equal-mass model LS 135135 shows only the m — 
and m — 2 modes dominating, with the addition of the "2-0" and 

1 These combination frequencies are also called bilinear coupling (or in- 
termodulation) components. 
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"2+0" combination frequencies. As in the previous models, there is 
again an /_ vs. "2-0" coincidence. 

For the strange matter EOS MIT60, we examined the equal- 
mass case MIT60 1111, where each individual component of the 
binary system has a mass of only 1.1 Mq (Fig. [7}. The reason is 
that for an individual mass of 1.35 Mq, an equal-mass binary sys- 
tem forms a black hole soon after merging. Notice that the density 
profile of the MIT60 models differs drastically from the case of 
the hadronic EOSs. It is very flat and terminates at a high value 
at the surface. Especially for the low-mass model of 1 . 1 Mq , the 
density profile is roughly uniform. This causes the oscillation prop- 
erties of the MIT60 models to differ considerably from those of the 
hadronic models. The GW spectrum of the MIT60 1111 model is 
still qualitatively similar to the previous hadronic models and one 
can still identify a triplet of frequencies /_ , fa and /+ . However, in 
this case, the frequency of the quasiradial mode is higher than for 
the hadronic models and, in fact, coincides with the frequency of 
the "2-0" nonlinear component. The latter differs, for this model 
only, from the /_ peak in the GW spectrum. It is possible that 
the /_, fa and /+ triplet is caused by the nonlinear interaction 
of the fa mode with a mode other than the quasiradial mode or 
that these are combination frequencies of higher order. Notice also 
that, for this model, the frequency of the m = 2 mode is twice 
the frequency of the quasiradial mode. Such a coincidence, which 
leads to a resonance between the two modes and to enhanced GW 
emission has also been observed in simulations of phase-transition- 
induced instabilities in rotating compact stars ( Abdik amalov et al.l 
l2009l : lDimmelmeier et al.l2009h . 

For the unequal mass MIT60 12135 model (Fig. [Tj, the fre- 
quency of the quasiradial mode is smaller than for the MIT60 1111 
model, while the frequency of the m = 2 mode is higher. As a con- 
sequence, the frequency of the "2-0" nonlinear component is not 
far from the /_ peak in the GW spectrum, although the agreement 
is not as good as in the case of the hadronic EOSs. One should 
keep in mind, however, that the post-merger remnant is evolving in 
time, contracting its radius and reaching higher densities. This is 
causing the different oscillation frequencies to change continually 
in time. One cannot, therefore, expect a perfect match between the 
measured peaks in the Fourier transform of evolved variables and 
the nonlinear features of the GW spectrum. 

We thus find that the GW spectrum of the post-merger phase is 
characterized by a triplet of frequencies that clearly coincides with 
the frequencies of the m — 2 oscillation mode of the fluid and its 
interaction with the quasiradial mode in all cases with mass in the 
range 1 .2 Mq to 1 .35 Mq . Only in the low-mass MIT60 1111 case 
the triplet has a different origin. 



remnant, could be used to distinguish a soft EOS from a stiff EOS. 
IShibata & Urvul J2002h also suggested that the post-merger GW 
spectrum is produced by non-axisymmetric oscillation modes of 
the merged object and arrived at empirical relations (for polytropic 
EOSs) for the two frequencies /_ and fa, without identifying their 
origin. Based on the dependence of these two frequencies on the 
parameters of the polytropic EOS, they suggested that these two 
freque ncies could be used t o constrain the stiffness of the EOS. 

In Shi bata et al.l 12005J) it was suggested that quasiradial oscil- 
lations modulate the post-merger signal and it was noticed that the 
difference in frequency between the two main peaks (i.e. the peaks 
we call /_ and fa here) is approximately e qual to the frequency 
of the quasir adial oscillation. Fur thermore, in lShibata & Taniguchil 
(2006) and in iKiuchi etalj|2009l) two strong sideband peaks (form- 
ing a triplet with the fa peak) were pointed out to exist in some 
cases and were associated with a modulation by large-amplitude 
quasiradial oscillations. The results of our mode analysis are in line 
with the above o bservations and e xpectations. 

In contrast, Bai otti et al. I d2008l) described the repeated appear- 
ance of a double core structure in a particular high-mass merged 
object as a dynamical bar-mode instability which develops and is 
quickly suppressed again several times. The double core structure 
appeared at a regular interval of 2 ms. According to our current 
results, this should be simply the period of the quasiradial oscilla- 
tion and instead of a dynamical instability the modulation is due 
to the "2- 0" nonlinear comp onent (in other, low-mass models pre- 
sented in lBaiotti et al .1 d2008l) . a bar persists for a long time and a 
possible association to a dynamical instability is worth further in- 
vestigation). 



5 COMPARISON TO PREVIOUS WORK 

lAllenetalJ ( ll999l) presented an initial study of binary neutron star 
mergers as a linear perturbation problem of an isolated star (a so- 
called close-limit approximation, in analogy to the corresponding 
approximation in binary black hole mergers). They showed that 
several fluid modes would be excited, but their study, being lin- 
ear, did not take into account nonlinear combination frequencies, 
as we do here in our fully nonlinear approach. 

In lZhuge et alj i 19941) the frequency peak /_ in the GW spec- 
trum was associated with the rotating barlike structure formed im- 
mediately after merging, while it was suggested t hat t he main fa 
peak is due to a low-order p-mode. Oechslin et al. (2002) proposed 
that the fa peak, being the quadrupole frequency of the post-merger 
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Figure 3. Model Shen 12135 - Top panel: Time evolution of 
the GW amplitude h+. Middle panel: Total (black), pre-merger 
(red) and post-merger (green) scaled power spectral density, com- 
pared to the Advanced LIGO and ET unity S N R sensitivit y curv es 
Harry & the LIGO Scientific Collaboration! bold ): Iffild et alj feOld) ). 
The distance to the source is assumed to be 100 Mpc. Bottom panel: Am- 
plitude of FFT for the time evolution of the pressure, p, in the equatorial 
plane. Several oscillation modes, as well as nonlinear combination frequen- 
cies (blue labels) are identified. 




Figure 4. Same as Fig. [3] but for model Shen 135135. 
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6 CONCLUSIONS 

We have studied the excitation of nonaxisymmetric oscillations 
in the post-merger phase of binary compact object mergers. Our 
analysis is based on general-relativistic simulations, using SPH 
for the evolution of matter, and we used a set of equal-mass 
and unequal-mass models, described by two nonzero-temperature 
hadronic EOSs and by one strange star EOS. We studied the os- 
cillations through Fourier transforms of the evolved matter vari- 
ables and identified a number of oscillation modes, as well as sev- 
eral nonlinear components (combination frequencies). The domi- 
nant m — 2 mode forms a triplet with two nonlinear components 
that are the result of a coupling to the quasiradial mode. A corre- 
sponding triplet of frequencies was identified in the GW spectrum, 
when the individual masses of the compact objects are in the most 
likely range of 1.2 to 1.35 M©. A specific frequency peak in the 
GW spectrum can thus be associated with the nonlinear compo- 
nent resulting from the difference between the m — 2 mode and 
the quasiradial mode. This association is especially strong in the 
case of hadronic EOSs and could be exploited in the case of fu- 
ture detections, in order to characterize the properties of the post- 
merger remnant. For this, it will be necessary to obtain accurate 
frequencies of the quasiradial and quadrupole oscillation modes for 
a large sample of theoretically possible post-merger remnants and 
construct empirical relations, depending on a few gross properties, 
such as the mass and radius of the star. Given the determination of 
two frequencies and of the total mass of the system from the inspiral 
signal, these empirical relations could be inverted to yield crucial 
information on the properties of high-density relativistic objects. 

It would further be interesting to test the sensitivity of our 
results to the spatial conformal flatness approximation employed 
here, as well as to assess the influence of magnetic fields on the 
gravitational wave spectrum. The latter effect can be expected to 
be small for realistic magnetic fields, although there could be a re- 
gion in the mass vs. magnetic field parameter spa ce in which the ef- 
fect co uld be measurable (for recent results, see iGiacomazzo et alj 
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